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_ ABSTRACT 

u 

• , Aims. We present the final public data release of the VLT/ISAAC near-infrared imaging survey in the GOODS-South field. The 

. survey covers an area of 172.5, 159.6 and 173.1 arcmin^ in the J, H, and A", bands, respectively. For point sources total limiting 

Ph' magnitudes of J = 25.0, H = 24.5, and K, = 24.4 (5cr, AB) are reached within 75% of the survey area. Thus these observations are 

significantly deeper than the previous EIS Deep Public Survey which covers the same region. The image quality is characterized by a 
, point spread function ranging between 0.34" and 0.65" FWHM. The images are registered to a common astrometric grid defined by 

the GSC 2 with an accuracy of ~ 0.06" RMS over the whole field. The overall photometric accuracy, including all systematic effects, 
adds up to 0.05 mag. The data are publicly available from the ESO science archive facility. 

Methods. We describe the data reduction, the calibration, and the quality control process. The final data set is characterized in terms 
of astrometric and photometric properties, including the PSF and the curve of growth. We establish an empirical model for the sky 
background noise in order to quantify the variation of limiting depth and statistical photometric errors over the survey area. We define 
^ ' a catalog of /fs-selected sources which contains JHKs photometry for 7079 objects. Differential aperture corrections were applied to 

\^ ' the color measurements in order to avoid possible biases as a result of the variation of the PSF. We briefly discuss the resulting color 

, distributions in the context of available redshift data. Furthermore, we estimate the completeness fraction and relative contamination 

fT^ . due to spurious detections for source catalogs extracted from the survey data. For this purpose, an empirical study based on a deep 

image of the Hubble Ultra Deep Field is combined with extensive image simulations. 

Results. With respect to previous deep near-infrared surveys, the surface density of faint galaxies has been established with unprece- 
dented accuracy by virtue of the unique combination of depth and area of this survey. We derived galaxy number counts over eight 
magnitudes in flux up to 7 = 25.25, H = 25.0, = 25.25 (in the AB system). Very similar faint-end logarithmic slopes between 
0.24 and 0.27 mag"' were measured in the three bands. We found no evidence for a significant change in the slope of the logarithmic 
galaxy number counts at the faint end. 



Key words. Cosmology: observations - Cosmology: large scale structure of the universe ~ Galaxies: evolution - Infrared: galaxies - 
Surveys 

1. Introduction facilities in space and on the ground, i. e. the NASA Great 

Observatories, the Spitzer Space Telescope, Hubble (HST), and 
The Great Observatories Origins Deep Survey (GOODS) aims chandra, ESA's XMM-Newton, the 8.2 m-aperture Very Large 
at combining the best and deepest data from X-ray through Telescopes at the ESO Paranal Observatory, the twin 10-meter 
radio wavelengths and making them publicly available to the ^^^^ Telescopes, the 8.2-m Subaru Telescope, and the 8-m 
community. The combined data enable studies of the distant telescopes of the Ge mini Observatory . For an overview of the 
universe in terms of the formation and evolution of galax- qqq^^ -^^^^ Dicki nson et al.l (l200l . iGiavalisco et al.1 
les and active galactic nuclei, the distribution of dark and lu- ^2004), and .Giavahsco (20051). 
minous matter at high redshift and the origin of extragalac- 

tic background radiation. Observations are centered on the two _ ,. 
target fields: the Hubble Deep Field North and the Chandra ,,,J° ^^udy the evolution of galaxies selection at near-infrared 
Deep Field South (CDF-S), each of them covering an area of ^IR) wavelengths has a number of advantages over optical se- 
160 arcminl The data are obtained using the most powerful ^^'^^1°"' *^ most significant one being t he smaH dependen ce on 

type out to redshitts ot about three (e. g. lCowie et alJll994 ). The 



Send offprint requests to: J. Retzlafi^, e-mail: jretzlaf@eso . org discovery of a population of distant galaxies based on NIR color 
* Based on observations using the Very Large Telescope at the selection (see, e. g.,|Franx et al.||2003|; |Daddi et al.||2004|) which 
Paranal Observatory of the European Organisation for Astronomical would have been missed by the Lyman break technique demon- 
Research in the Southern Hemisphere, Chile, under ESO programmes strates the importance of deep NIR imaging for the understand- 
168.A-0485, 64.O-0643, 66.A-0572, 68.A-0544, and 73.A-0764. ing of the history of mass assembly in the universe. 
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While deep ground-based optical imaging has become 
largely obsolete due to the GOODS ACS data set, high-quality 
NIR observations had to be executed from a ground-based facil- 
ity because the small field of view of NICMOS on HST would 
have made NIR imaging of the GOODS region overly costly. 
However, NIR imaging remains a key requirement for the scien- 
tific goals of the GOODS program because sampling of the NIR 
spectral range is crucial if spectral energy distributions have to 
be fitted in order to derive accurate photometric redshifts and 
stellar population parameters of galaxies. ESO has dedicated a 
substantial amount of observing time in order to support the 
community with imaging and also VLT spectroscopy for the tar- 
get field CDF-S (Renzi ni et al. 2003.) . 

There were four pubUc incremental releases of the GOODS 
ISAAC data: version 0.5 (April 2002), version 1.0 (April 2004), 
version 1.5 (September 2005), and version 2.0 (September 
2007). In this publication we present the final version (2.0) of the 
fully reduced and calibrated images which have been obtained 
with the Infrared Spectrometer And Array Camera (ISAAC) in 
the J, H, and K^^ bands. This data set was released via the ESO 
science archive facility as part of the ESO/GOODS projectQ 

The paper is organized as follows. Sect. |2] outlines the ob- 
servations. Sect. [3] describes the image data reduction and cali- 
bration process. The final survey images, their properties and an 
assessment of the photometric calibration are given in Sect.H] In 
Sect. |5] we describe how photometric source catalogs were de- 
fined and discuss their characteristics in term of color-magnitude 
and color-color diagrams. In Sect. |6] first, catalog completeness 
and contamination is discussed, and thereafter deep galaxy num- 
ber counts in the three survey bands are presented. Finally, we 
conclude in Sect. [7] 

Note that throughout this work magnitudes ar e expressed in 
the AB photometric system dOke & Gunnlll983h unless other- 
wise noted. 
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Fig.l. TiHng of the CDF-S region with 2.5' x 2.5' ISAAC 
fields displayed over the HST/ ACS GOODS z-band mosaic 
(iGiavalisco et al.l2004l) . The fields on the very top, FOl and F02, 
are covered in J and K^, the other 24 fields (F03 through F31) in 
all three bands J, H, and K^^. The dashed diamond centered on 
field F14 is the footprint of the deep UDF/ISAAC Ks image. 



2. Observations 

The Infrared Spectrometer and Array Camera (ISAAC) mounted 
on the first of the four 8.2-m VLT unit telescopes at Cerro 
Paranal is equipped with a 1 024 x 1 024 pixels HgCdTe Rockwell 
Hawaii array having a pi xel scale of 0.148", co rresponding to a 
2.5' X 2.5' field of view ( Moorwood et alJ[T99 9). Fig. [1] shows 
how the ISAAC pointings were laid out to assemble a mosaic of 
the CDF-S region. A contiguous area of 24 fields, 150 arcmin^ 
respectively, are covered in the J, H, and Ks bands, plus 2 addi- 
tional fields (FOl and F02) at the top of the survey area which 
have J and but no H data. In the early stages of the project, 
in an attempt to maximize the overall survey efficiency, fields 
F25 and 26 were re-arranged and replaced with F25n and 26n. 
The Ks band data which was obtained for fields F25 and 26 has 
been included in the survey data processing in the same way as 
the other fields despite their shallowness. Later on, this auxiliary 
data proved to be very valuable for the purpose of validation of 
the internal photometric consistency (Sect. [43] ). See Table[T]for 
the field coordinates. 

The GOODS ESO ISAAC program was originally planned 
to reach 5a" limiting magnitudes of 25.2, 24.7, and 24.4 in J, 
H, and K^, respectively, in order to roughly match pre-launch 
expectations for Spitzer IRAC sensitivity at 3.6 to 8 jum. The 



' To be downloaded from 
http : //archive . esc . org/cms/eso-data/data-packages/ 
goods-isaac-final-data-release-version-2-0 



intent was to provide NIR flux measurements for the large ma- 
jority of IRAC-detected galaxies, as well as images with higher 
angular resolution in order to facilitate the deblending of sources 
in crowded regions of the lower-resolution IRAC data. In prac- 
tice, Spitzer IRAC in-flight performance at 3.6 and 4.5 /im 
significantly exceeded the pre-launch expectations, making the 
GOODS data at those wavelengths substantial deeper than any 
ground-based NIR program could reasonably match. Given the 
instrumental sensitivity, a total integration time per tile of 3.5, 
5, and 6 hours was estimated to reach the intended depth. To 
achieve accurate sky background subtraction which is most cru- 
cial for deep NIR imaging, the commonly practiced jitter imag- 
ing technique was used. A jitter box size of 25" was cho- 
sen within which the control system automatically off'sets the 
telescope between subsequent integrations. Detector integration 
times (DITs) between 10 and 30 s were used, according to the 
typical sky brightness (darker at J, brighter at H and Ks). The 
number of detector integrations (NDIT) was chosen so that to- 
tal integration times (DIT x NDIT) between 60 and 180 s were 
used. The survey was executed in observation blocks (OBs) each 
of which comprising the integration of a single field in one filter 
for a total integration time of normally 3600 s or 1800 s, and, 
thus, resulting in between 15 and 40 images. The total actual in- 
tegration time of the data that were combined in the final survey 
images amounts to 360 hours. Given that some data were dis- 
carded and taking into account observing overheads the whole 
program was allocated about 500 hours of observing time. The 
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Table 1. Nominal center positions of the GOODS/ISAAC survey 
fields. 
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observations have been executed in service mode in 216 nights 
between October 1999 and January 2007. 

The data were obtained under the ESO large programme 
168.A-0485, led by C. Cesarsky, in direct support to the GOODS 
project. Data covering four ISAAC fields in J and bands were 
also drawn from the ESO programmes 64.O-0643, 66.A-05 72 
and 68.A-0544, led by E. Giallongo (see lSaracco et alj|200ll) . 
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Final imaging 
data products 



Fig. 2. Schematic flow of the data reduction and calibration pro- 
cess from raw data to final, fully calibrated survey data products 
(from top to bottom). Processing steps are represented by rect- 
angles, input/output data by parallelograms, and thick arrows in- 
dicate the image data flow. Numbers represent the cardinality of 
each data set. Refer to the text for a detailed description. 



3. Image data reduction and calibration 

3A. Overwew 

The overall survey data reduction and calibration process is 
schematically illustrated in Fig.|2l To begin with, the ESO/MVM 
data reduction system was used to process the raw jitter mode 
observations using the required calibration files to correct for in- 
strumental signatures and to produce OB images, that is one co- 
added image with associated weight map per OB being astromet- 
rically registered to the predefined astrometric grid of the survey. 
The fully automated processing of the entire data set of 13964 
raw jitter mode frames and 20699 calibration files resulted in 
426 astrometrically calibrated OB images (Sect. [372)13. 31 ). Then, 
the OB images were characterized in terms of PSF width, pho- 
tometric uniformity, and depth in order to identify and remove 
images of low-quality which would have degraded the quaUty of 
the final survey products (Sect. [341 ). 

For the photometric calibration of the OB images, at first, 
the raw standard star observations were reduced to instrumen- 
tal photometric zero points using the master flatfield calibrations 
generated by ESO/MVM. From this data set, a subset of "best" 
zero points was constructed. Then, photometric scalings of OB 



images belonging to the same tile and between adjacent tiles 
were measured and combined with the zero points to determine 
- in a robust way - an accurate photometric calibration for all 
OB images equalized over the whole survey field (Sect. 13.5b . 
That followed, the photometrically calibrated OB images were 
co-added to produce the survey tiles which subsequently under- 
went a final step of image quality characterization. 

3.2. Data processing 

The ESO/MVM software, version 1.3.4, has been used for im- 
age data reduction. ESO/MVM follows well established proce- 
dures for the reduction of NIR imaging observations taken in 
jitter mode by featuring a two-pass scheme to mask out sources 
prior to the final background estimation. In the following, we 
will outline the main reduction steps executed by ESO/MVM. 
Comprehensive documentation of ESO/MVM including all the 
reduction steps and algorithmic details being impleme nted is be- 
yond t he scope of this publication and may be found in lVandam3 
(12004 which serves as reference documentation and describes in 
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length the implemented algorithms, parameter settings, instru- 
ment specific configurations, and sample application^ 

First, the raw calibration data, namely dark frames and twi- 
light sky flats, were combined into master calibration files and 
used to correct the raw scientific images for the basic instru- 
mental signatures. No attempt was made to correct for detector 
nonlinearity. However, the overall effect is expected to be rather 
small because the scientific OBs for the survey program and, 
by the same token, the OBs for photometric standards had been 
generally designed with the objective to operate the detector in 
the linear regime under nominal conditions. The magnitude of 
possible nonlinearity effects on the final results is discussed in 
Sect. [3331 

In a few cases, overly long OBs were split in two blocks, or, 
particularly short, partially completed, OBs were merged prior to 
data reduction in order to optimize the quality of the reduction. 
Sky background images were computed from groups of between 
13 and 15 consecutive jitter images (depending on filter) us- 
ing sigma-clipped pixel-by-pixel image combination. Each sci- 
ence image was sky-subtracted using the linear interpolation in 
time of the two corresponding successive sky background im- 
ages. Possible transparency variations from image to image were 
monitored by means of the signal-to-noise ratio (SNR) based on 
which outliers with exceptionally low SNR were automatically 
discarded. An individual rescaling of the images has not been 
done. 

Then, for each OB, the individual background-subtracted 
images were astrometrically registered to each other with sub- 
pixel accuracy and co-added (see below) to form a preliminary 
version of the OB images. The purpose of these first-pass im- 
ages is the creation of masks that mask the astronomical sources 
in order to improve the sky background computation in the sec- 
ond pass. To this end all sources exceeding ~ 3 times the local 
RMS noise were detected on the PSF-convolved image and - to 
also exclude the wings of the source profile - each source's re- 
gion was artificially enlarged by a factor of 2 (linearly) before 
creating its mask. Thus, in the second pass, the sky background 
images were re-computed exclusively from image pixels that be- 
long to background regions whereas the rest of the procedure 
was repeated unaltered. 

3.3. Astrometric calibration 

The astrometric calibration is based on a dense reference catalog 
which was generated by the GOODS team from a deep R-band 
image of the CDF-S and was also used for the production of the 
GOODS/ACS image mosaics (Giavahsco et al. 2004). The im- 
age was obtained with the Wide Field Imager mounted at the 
2.2-m MPG/ESO telescope at La Silla, and astrometrically cali- 
brated using the Guide Star Catalog (GSC 2). Each OB image of 
the GOODS/ISAAC survey was astrometrically registered using 
the reference catalog whereas image distortions were modeled 
with a 3'^'' order polynomial. This resulted in an internal astro- 
metric accuracy between 0.05" and 0.06" RMS as measured for 
sources brighter than 20 mag in the final mosaic images in all 
three filters. The relative registration between bands is accurate 
to -0.03". 

For the astrometric grid of the OB images the same projec- 
tion as for the GOODS/ACS images has been adopted and a 
pixel size of 0.15" has been chosen so that one GOODS/ISAAC 



^ The full software package may be downloaded from 
http : / /archive . eso . org/cms/eso-data/ data-packages/ 
eso-mvm-software-package 



pixel subtends exactly a block of 5 x 5 GOODS/ACS pixels. 
The individual jitter images were resampled to this grid using 
the Lanczos-3 interpolation kernel. Pixel-to-pixel noise corre- 
lation, the interpolation process has introduced, is quantified 
in Sect. 13.41 During the process of image co-addition, bad- 
pixel masks were taken into account and each contribution was 
recorded pixel-by-pixel to build up respective weight maps. 

The comparison of the final mosaics with the calibrated data 
from the Hubble Space Telescope (HST) Advanced Camera for 
Surveys (ACS) in its version 1 . 1 incarnation yields astrometric 
offsets of typically -0.1" RMS across the entire area which can 
be attributed primarily to zonal residuals that are known to be 
present in the astrometric calibration of the HST/ACS mosaics. 

3.4. Quality control of OB images 

For each OB image the quality was quantified in terms of the 
FWHM of stellar sources ("seeing"), the sky noise as a function 
of aperture size, and the photometric scaling of images belong- 
ing to the same tile and filter Then, based on the inspection of 
these quality parameters with respect to the whole sample, im- 
ages of significantly low quality were identified and sorted out, 
followed by a visual inspection, to ensure that no corrupted im- 
age slipped through. 

In order to determine the average FWHM of stellar sources, 
we made use of a list of ca. 400 sources that appear point-like 
in the GOODS/ACS z-band image. Based on this list we pre- 
selected appropriate candidate sources for the PSF characteriza- 
tion. 

The sky noise was measured for each OB image in a similar 
fashion as described bv lLabbe et al.l (l2003b . that is, first, a noise- 
equalized image was created by multiplication with the square 
root of the weight map. Then, the image was segmented into 
background and sources using a ~ 2.6cr significance threshold 
plus a 5 pixel wide safety margin around each source segment. 
Then, the background flux was sampled with non-overlapping 
circular apertures with diameters between 0.5" and 4.0" and a 
linear relation of the form {a + bs)s was fitted to the dispersion to 
obtain the two parameters a and b. Here, the linear aperture size, 
s, is defined as the square root of the aperture area in pixel units. 
It has been verified that the adopted parameterization describes 
the actual data adequately. Therewith, the sky noise, that is the 
statistical fluctuation of the flux measured within an aperture of 
size s at any point of the image can be expressed by the Gaussian 
dispersion 

Co 

oskyis) = — =(a + bs)s, (1) 

where ctq is the RMS of all background image pixels, and wj is 
the relative weight of pixel j. This empirical model of the sky 
noise accounts for the pixel-to-pixel correlation which results 
from the preceding image interpolation step and which breaks 
the simple scaling with aperture scale cTsj^y oc s. According to 
the interpolation kernel, the correlation length is not more than 
a few pixels. To quantify the amount of this small-scale con- 
tribution, the total a + b is typically about 1.15, that is, at the 

^ For the relative weight, Wj, a normalization to unity at the median 
pixel has been adopted, i. e., w, = lVy/median(lV,) with Wj denoting the 
inverse variance weight map and the median being taken of all pixels of 
the given image having Wj > 0. Thus, cro refers to the median pixel. 
This notion, using the index zero, has been adopted throughout the text 
for quantities scaling with the weight map, namely the limiting depth 
and the gain conversion factor. 
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Fig. 3. Variation of the photometric zero point of the 
GOODS/ISAAC OB images from October 1999 until January 
2007. Filled symbols correspond to OBs with flux standard star 
calibrations, open symbols refer to OBs whose zero points have 
been inferred by photometric scaling. Zero points have been con- 
verted to unit airmass to produce this plot. However, this correc- 
tion is negligibly small compared to the gross variation. 



scale of a single pixel, i. e. s = 1, the direct estimation ctq un- 
derestimates the flux dispersion by about 15%. The other con- 
tribution which affects slightly larger scales is due to imperfect 
background subtraction, that is residual spatial modulations of 
the image background. If the night sky background is highly 
variable with a time scale being shorter than the one implicit 
to the averaging process that is part of the image data reduction, 
then such residuals may remain - mostly apparent in the form of 
ripples or speckles. Not surprisingly, the parameters a and b are 
sensitive indicators for the presence of such residuals. Therefore, 
we have employed them to identify and reject OB images which 
were manifestly aff'ected above average. 

In the course of the quality control process, 15 OB images 
were rejected because of the seeing being worse than 0.85" 
FWHM, nine for being very noisy and showing particularly 
strong residuals in the sky background {b > 0.15), five for be- 
ing exceptionally "shallow" with the photometric zero point of 
more than 0.3 mag below average, and three for other reasons 
such as corrupted data. In total, 394 OB images, which account 
for 93.7% of the total integration time of reduced OBs, passed 
the quality control criteria and underwent the following steps of 
photometric calibration, image stacking, and final characteriza- 
tion which we will describe in the next section. 



3.5. Photometric calibration 
3.5.1. Photometric zero points 

As the GOODS/ISAAC program was spread over a large time 
frame, the substantial variation of the instrumental photometric 
zero point on various time scales is not surprising (Fig. [3]). These 
variations are caused by instrumental interventions, the steady 
degradation of the mirrors' reflectivity, mirror re-coating events, 
and, finally, the variation of atmospheric transparency. The sig- 
nificantly low instrumental zero point at the beginning of VLT 
operations, for instance, when the J and Ks data for tile F16 were 
being obtained, is likely due to accumulated dust on the mirror of 



UTl caused by the ongoing construction of the other UTs at that 
time. In order to establish an accurate and consistent photomet- 
ric calibration across the whole survey area, we have combined 
photometric zero points (ZP) obtained from standard stars with 
relative photometric scalings between survey tile images. 

Instrumental ZPs were obtained using faint NIR standard 
stars which are being observed on a nightly basis as part of 
the observatory's instrument calibration p lan. Most of the flux 
standa rds were drawn from Persson et al.Tfl998l) and lHunt et al.l 
()1998'), complemented by a small number of UKIRT standards 
CHawarden et al. 2001). 

The available data of photometric flux standard observations 
was carefully selected to minimize the potential bias due to de- 
tector nonlinearity. For each image of a photometric standard 
star the peak signal including the background was used as qual- 
ity indicator and outliers exceeding 1 5000 ADU which corre- 
sponds to a formal nonlinearity of 1 % dAmico et alj|2002h were 
eliminated. The remaining set of photometric calibrations from 
which the ZPs were finally determined is dominated to 85% by 
data with peak fluxes below 10000 ADU which corresponds to a 
formal nonlinearity of ca. 0.5% while the remaining 15% of the 
data have peak fluxes above 10000 ADU. Additionally, we have 
examined more than 5000 individual photometric calibrations in 
J, H, and band as a function of the maximum signal level to 
empirically check for a possible systematic trend in the derived 
ZPs. Up to a signal level of 15000 ADU, we found no evidence 
that ZPs are systematically biased low. Merely beyond 20000 
ADU, i. e. for data that was finally not used for calibration, a de- 
cline by ^0.03 mag becomes apparent in all three bands. Hence, 
we can conclude that the ZPs which were eventually used to an- 
chor the survey's photometry are affected by detector nonlinear- 
ity by less than 0.01 mag even in the worst case. 

We started out by selecting all the standard star observations 
within an interval of plus or minus 2 nights around each science 
observation and derived a total of 776 zero points - by far more 
than what went into the final photometric solution (see below). 
ZPs were computed from the instrumental flux me asured within 
an aperture of 10" diameter in compliance with iPersson et al.) 
(1998). 

Since flux standards had not always been observed imme- 
diately before or after each science integration, we adopted the 
following scheme for the association of ZPs aiming at a practi- 
cal compromise between conservatively minimizing the effect of 
possible transparency variations and ending up with a sufficient 
number of ZPs. For each OB the ZP whose flux standard had 
been acquired closest in time to the science observation within a 
time interval of less than 2 hours and an airmass difference be- 
low 0.5 was selected. If these conditions were not met by any 
ZP, the respective OB was not assigned a ZP but the photomet- 
ric calibration was purely inferred from the scaling relative to 
other OBs belonging to the same or to neighboring tiles (see be- 
low). Using these criteria, 210 out of the 394 OB images were 
associated with 165 ZPs, 50% of which were obtained in a time 
interval of plus or minus 45 minutes with respect to the corre- 
sponding science observation. 

Then, we have inspected each associated ZP in the context 
of all the other ZPs obtained within an interval of 5 nights so 
as to identify non-photometric conditions and, hence, we have 
dismissed 10 ZPs for being low by S 0.05 mag. 

The resulting 155 instrumental ZPs were converted into 
ZPs for 199 OB stacks using the atmospheric extinction coef- 
ficients for the ISAAC instrument, 0.09, 0.04, and 0.06 mag per 
unit airmass for J, H, and K^, respectively, as determined by 
iMason et al.l (l2008l) . 
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Fig. 4. Residuals of photometric zero points of OB images as de- 
rived from flux standards with respect to the final global photo- 
metric solution for J, H, and (from top to bottom). For each 
passband, the total number of OB images having photometric 
zero points based on flux standards, A^, and the RMS of the resid- 
uals about the mean, cr, is reported. Gaussian distributions with 
the same mean and dispersion as the actual data are superposed. 



The match between the ISAAC filters (J, H, and K^) and 
those used to estab lish the faint IR standard star system of 
iPersson et alj (Il998l) is quite good and one expects color t erms 
which differ from by less than 0.01 (lAmic o et aLll2002 l). As 
the color transformation between ISAAC magnitudes and those 
of LCO have never been experimentally verified, and our data 
does not allow for it either, we have chosen not to apply any 
color correction. 

3.5.2. Global photometric calibration 

To determine a survey-wide photometric solution per filter, pho- 
tometric scalings between individual OB images were com- 
puted. As it was straightforward to determine the relative flux 
scaling for images belonging to the same tile and filter with sub 
per-cent accuracy once PSF matching had been done, the scal- 
ing for neighboring tiles turned out to be difficult to measure, 
because, usually, the overlapping area was too small to find a 
sufficient number of high signal-to-noise sources for accurate 
photometry. Therefore, we have employed public data from the 
ESO science archive which covers the survey area in the same 
bands in order to establish accurate relative photometric calibra- 
tions between adjacent surve y fields. The data had b een obtained 
using the SOFI instrument (IMoorwood et al.iri998l) mounted on 
the New Technology Telescopy at La Silla. SOFI's field of view 



of about 4.9' x4.9' guarantees sufficient overlap between ISAAC 
and SOFI images in general, while other instrument characteris- 
tics are similar by construction, in particular, the instrumental 
response. The J and observations that were use d belong to 
the in frared part of the ESO Imaging Survey (e. g. lOlsen et al.l 
I2006h . the H band data were orig inally obtained w ithin a pro- 
gram led by D. Rigopoulou (see M ov et al] l2003h Pl The lay- 
out of the JHKs SOFI observations in the CDF-S generally 
includes some overlap of adjacent images, and, additionally, a 
series of low-exposure "calibration images" arranged in a pat- 
tern being offset to the deep integrations had been acquired in 
J and K^, so that accurate photometric scalings could be mea- 
sured not only between ISAAC and SOFI but also between adja- 
cent SOFI images. To this end, we have processed the raw SOFI 
data in the same way as we did with the ISAAC data, employ- 
ing ESO/MVM for the data reduction followed by careful ex- 
amination of the image quality (Sect. 13.41 ). After all, 39 best- 
quality images in terms of seeing, noise, and background homo- 
geneity were used in J, H, and Ks, plus 19 shallower "calibra- 
tion images" in J and Ks- To measure the photometric scaling 
between two images, at first, the image with the better seeing 
was smoothed with a Gaussian kernel to match the PSF of the 
other image with poorer seeing. Because typically the ISAAC 
seeing is better than the SOFI seeing (0.5" vs. 0.7" median), 
for scalings between ISAAC and SOFI images in the major- 
ity of cases the ISAAC image is smoothed to match the PSF 
of the SOFI image. Then, instrumental magnitudes of all iso- 
lated, high signal-to-noise (SNR>20) sources were measured 
using SExtra ctor's auto-scaling apert ure magnitudes in double- 
frame mode ( Bertin & Arnout si ll 9961) and the magnitude differ- 
ences were averaged to obtain the relative photometric scaling. 
In this way, photometric scalings with respect to SOFI images 
could be measured with a typical uncertainty between 1 and 2%, 
slightly depending on the filter We have minimized a possi- 
ble bias due to the relative nonlinearity of the two instruments 
as far as possible. Objects being relatively bright for ISAAC 
were discarded based on their peak flux corresponding to an 
effective flux cut at ~ 16.5 mag (AB). Furthermore, the pho- 
tometric scaling was computed by averaging over all suitable 
sources, i. e. typically more than 10 sources, stars and galax- 
ies, between ca. 18 and 20 mag (AB) contribute to the resulting 
scaling instead of being based on a few bright stars only. For 
objects fainter than 17 mag, ISAAC'S detector nonlinearity is 
not an issue (Sect. [333T l. As another precaution, we visually in- 
spected the differential ISAAC-SOFI photometry but did not see 
any case of flux dependent bias. Therefore, we can exclude that 
the measured photometric scalings between ISAAC and SOFI 
data are biased significantly, that is by more than 0.01 mag. That 
followed, we have a pplied a global ;(f^-minimization technique 
dKoranvi et al.ll 199*81) . to find, for each passband, tiie photomet- 
ric solution adjusted across the whole survey. 

The residual photometric differences of the photometric ZPs 
as given by the flux standards with respect to the global pho- 
tometric solution were examined to identify implausible ZPs. 
Consequently, 5 ZPs were removed, and the final photometric 
solution was calculated based on 155 photometric ZPs and 501 
photometric scalings. This means that on average, more than one 
photometric ZP contributed to the final photometric solution per 
OB image, namely 1.4, 2.5, and 3.2 for J, H, and Ks, respec- 
tively, thereby reducing the impact of possible errors of individ- 
ual photometric ZP measurements by virtue of averaging. As an 



The observations have been obtained under ESO programmes 
60.A-9005, 164.O-0561, 66.A-0451, and 68.A-0375. 
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Fig. 5. Reproducibility of flux measurements on OB images as a 
function of AB magnitude. Each data point corresponds to one of 
56 individually selected test sources measured in one of the three 
filters, J, H, and K^, discriminated by different plot symbols, 
square, diamond and circle, respectively. The root-mean-square 
(RMS) flux deviation from the mean was computed from 3 to 
30 measurements on the photometrically calibrated OB images. 
Solid symbols denote unresolved sources, while open symbols 
correspond to extended sources. For comparison the Icr errors 
of the photometric solution of the OB images (Sect. 13.5.21 ) are 
shown as horizontal dotted lines. 



additional check, we have verified that the residuals are not cor- 
related with observational or instrumental parameters (seeing, 
sky background level, DIT). The photometric residuals appear 
to be normally distributed, and indicate an internal photometric 
accuracy of better than 0.02 mag RMS from tile to tile in all 
three bands (Fig. |4|. The individual inspection of a set of non- 
saturated, isolated stars has revealed in a few cases photometric 
discrepancies of 3-5% between SOFI and ISAAC images which 
were not attributable to photon noise but are presumably due to 
systematic errors left over after flat field correction. 

Finally, in order to convert from the Vega to the AB sys- 
tem, AB corrections of 0.9603, 1.426, and 1.895 mag for J, 

H, and K^, respectively, were applied. These numbers were 
obtained from the ESO Magnitude-to-flux-converte:3 (version 

I. 07) which makes use of the instrumental transmission curves 
of the ESO Exposure Time Calculators. 

3.5.3. Photometric uncertainties 

In order characterize the effect of detector nonlinearity on the 
final survey data we have performed photometric checks of in- 
dividual objects on the calibrated OB images as a function of 
source properties (flux, extent), observational conditions (see- 
ing, background brightness), and instrumental parameters (DIT). 
To this end five survey tiles were selected within which aperture 
fluxes (4" diameter) of 56 isolated objects were measured in J, 
H, and from between 3 and 30 OB images. Fig.|5]summarizes 
the results of this study as a function of source flux. The increas- 
ing RMS flux deviation towards the bright end indicates that de- 
tector nonlinearity starts to take eff'ect for unresolved sources at 
total magnitudes between 16 and 16.5 AB. The comparison with 



Accessible at http : // archive . eso . org/ apps/mag2f lux 



the calibrated SOFI images, which are expected to be signifi- 
cantly less affected by nonlinearity than the ISAAC data, allows 
to directly estimate the net flux bias for these objects. On aver- 
age the ISAAC photometry underestimates the flux by between 
0.02 mag at ~ 16.5 mag and ca. 0.05 mag at 14.7 mag AB irre- 
spective of the pass band. The maximum flux bias of 0.07 mag 
was found for the H flux of the brightest star in the test sample 
(H = 14.7). Note that there are actually very few objects that 
bright in the whole survey, namely 5 having < 15 and 23 
having < 16.5, almost all of them being stars with the excep- 
tion of two bright galaxies of - 16.1 and 16.5. For the bulk 
of objects having total fluxes of 17 mag AB and fainter, no indi- 
cation for nonlinearity effects has been found. In fact, the overall 
trend of the flux uncertainty, which resolved and unresolved ob- 
jects seem to follow in the same way, can be explained purely in 
terms of statistical errors (see also Sect. 14. 2t . 

Fig. |5] also indicates the errors of the photometric solution to 
illustrate the level of systematic photometric errors in proportion 
to the statistical errors. For instance, between ca. 16 and 18 mag, 
H and photometry is typically limited by the systematic un- 
certainties intrinsic to the photometric solution, or, photometry 
fainter than J ~ 20 becomes dominated by statistical errors. 
Remember that the data points in Fig. |5] refer to the Icr fluc- 
tuations of photometric measurements on OB images, whereas 
the final co-added survey images are expected to have reduced 
photometric errors, typically by a factor of between 2 and 2.5. 

Exceptionally high sky background brightness in combina- 
tion with a comparatively large detector integration time (DIT) 
has resulted in a few OBs in H and for which the detector was 
not operated in the linear regime. In fact, the most extreme cases 
are 8 out of 294 OBs in H and for which the sky background 
signal induces a nominal nonlinearity of between 2 and 3%. To 
compensate for this effect which is effectively resulting in a loss 
of sensitivity, or a lowered photometric ZP, the respective OBs 
were not directly anchored to ZPs but their photometric calibra- 
tion was inferred from relative photometric scalings with respect 
to the other OBs of the same tile and filter (see Sect. 13.5.2b . 
Finally, we have verified that the photometric measurements on 
the calibrated "nonlinear" OBs are indistinguishable from the 
other data. 

To estimate the overall error, we linearly add up the individ- 
ual contributions, that is the internal accuracy of the photometric 
solution (< 0.017 mag), the residual flat field error (0.025 mag, 
Icr), and the possible nonlinearity bias 0.01 mag) for objects 
not brighter than 17 mag (AB), resulting in the total photometric 
error of up to 0.05 mag (Icr). 

3.6. Final image co-addition 

The co-addition of the 378 OB images into 78 survey tile im- 
ages marks the final step of the image data reduction process (cf. 
Fig. 12]). There is no resampling involved in this step because the 
OB images are already defined on the final astrometric grid. 

At first, OB images were reseated to the fiducial ZP of 
26.0 mag (AB), which in fact is quite close to the typical value 
for all three bands (see Fig [3j. Relative weights between OB 
images contributing to the same tile were chosen to be propor- 
tional to the inverse variance of the flux integrated within the 
circular aperture that maximizes the total limiting magnitude for 
point sources (cf. Sect. 14.2b . The weight map associated with 
each OB image was employed to consistently accounted for the 
positional dependence of the variance. By comparison with dif- 
ferent weighting schemes, we have experimentally verified that 
the adopted scheme indeed optimizes the resulting image depth 
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Fig. 6. Limiting magnitude (5cr, total for point sources) of the final survey tiles as a function of (a) total integration time, and, (b) 
PSF size ("seeing"). Panel (c) shows the limiting magnitude corrected for the variation of the instrumental zero point and scaled by 
exposure time and seeing so that the dependencies formally cancel out (see text). The survey tiles are marked with square, diamond 
and circle symbols as indicated in the legend. The triangle symbol represents the UDF/ISAAC A's image. Dotted lines denote the 
median depth of the GOODS/ISAAC tiles. Histograms of seeing and limiting magnitude are marked by thick line (blue), moderately 
thick line (green), and normal line (red) plus grey shade for J, H, and K^, respectively. 



for point sources as intended. The associated final weight maps 
are inverse variance maps and were rescaled to the actual, empir- 
ically determined pixel-to-pixel variance of the sky background. 
For the final survey images the gain factor. Go, that is the num- 
ber of detector electrons per data unit, is reported in Table|2] last 
column. As Go refers to the median pixel, the effective gain at 
pixel j can be obtained by scaling with the relative (i. e. median- 
normalized) weight, Wj, according to Gj - GoWj. 

In addition to the individual image tiles, the data release also 
includes mosaics of the co-adjoined tiles as single FITS files in 
J, H, and bands, as well as the corresponding weight maps. 
The net area is 172.5, 159.6 and 173.1 arcmin^ in J, H, and Ks, 
respectively. The WCS information and accuracy of the individ- 
ual tiles is preserved in these mosaics. A uniform ZP of 26.0 mag 
can be used (e. g. with SExtractor) across the entire field, how- 
ever, it is important to note that the PSF varies from tile to tile 
within each mosaic. In the absence of proper aperture corrections 
or PSF matching procedures, this would lead to biases when cre- 
ating multi-color catalogs. A possible procedure for coping with 
the PSF variations without resorting to image convolution is out- 
lined in Sect. 15.21 



4. Final survey images 

4. 1 . General properties 

In the following, we summarize and discuss the properties of 
the final survey images as presented in Table |2] The recorded 
parameters are: total integration time, total number of raw im- 
ages which were combined to form the final product, the period 
during which observations where conducted, the FWHM of the 
PSF ("seeing") in arcseconds, the 5cr limiting magnitude correct 
to the total flux assuming a point source profile ("image depth"), 
the aperture diameter to which the depth refers, the 3 parameters 
of the noise model ctq, a, b, and the effective gain. 

Summing the actual exposure time from the data that were 
combined in the final stacks yields a total amount of integra- 



tion time of 359.9 hours. Regarding the distribution of integra- 
tion time that went into the final images, one notices many tiles 
standing out with respect to the nominal values (Fig.|6h). This is 
because OBs which were rejected in the course of the final qual- 
ity control process could not be re-scheduled for observations 
leading to reduced integration time for those tiles with respect 
to the nominal values for the survey. For instance, in the most 
extreme case (FQAKf^), three out of six OBs were rejected - two 
OBs for high noise, one for bad seeing. The two "supplemen- 
tary" tiles, F25 and F26, are an exception in the sense that they 
have just a fractional amount of observation time by design. On 
the other side, there are tiles with integration times in excess of 
the nominal values, most notably Fll and F16 in K^, which is 
primarily due to the fact that data from several observing pro- 
grams have been combined. 

The PSF of the final images, quantified by its FWHM using 
the same methodology as detailed in Sect. 13.41 varies between 
0.34" (F23//) and 0.65" (F15//) with a median value of 0.48" 
and appears to be independent of the pass band (Fig. |6j3). We 
have checked that the final FWHM does not show any correla- 
tion with the number of raw images or the total integration time, 
which indicates that the image reduction process does not de- 
grade the PSF quality. Generally, we have found a low level of 
PSF anisotropy. SExtractor-measured ellipticities are typically 
smaller then 0.05 - sometimes even clearly below 0.03. 

Based on the analysis of the curve of growth of unsaturated 
isolated stars using circular apertures between 0.7" and 10" di- 
ameter, we have quantified the aperture corrections for point 
source photometry, Snig^p, for each tile (Table |3ll. 10" has been 
adopted as the reference aperture for the practical reason, that 
this turns out to be the upper limit out to which the correction 
can be traced in case of most favorable conditions, i. e., if suf- 
ficiently bright stars are present in the image. In other cases - 
when appropriate stars are lacking - a power-law extrapolation 
of the curve of growth had to be applied to obtain the data for 
the largest apertures. For all images the curve of growth falls 
significantly below the 1 % level at the end which justifies to 
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consider the 10" flux as the total flux, and to infer the total 
magnitude, Wtot, from the aperture magnitude, m^p by means of 
mtot - ni^p - <5»Jap as is done throughout the rest of the work. 

4.2. Sky background and limiting magnitude 

First, we have evaluated the sky background noise using the 
same methodology as detailed in Sect. 13.41 That is, for each im- 
age, we have obtained a model for the sky noise according to 
Eq. [T]in terms of the 3 parameters ctq, a, and b. Below, we will 
also employ these numbers to estimate flux errors (Sect. 15. 11 1. 
Then, the total magnitude for point sources being associated 
with the sky background fluctuation was computed taking into 
account the aperture correction, (5map, and the photometric ZP, 

miin, = -2.5 log (Tsky - dm^p + ZP, (2) 

and its maximum with respect to the aperture diameter was iden- 
tified as the point source limiting magnitude of the given image 
at 1 cr significance level. As customary, we refer to the 1 .747 mag 
brighter 5cr limiting magnitude in the following. The variation 
across the image is given by the weight map, wj, according to 
'Wiim = + 1 -25 log Wj, whereas m{||^ refers to the median pixel 
by convention. The median limiting magnitude and the diameter 
of the corresponding circular aperture are listed in Table |2] 

The primary dependencies of the limiting magnitude of the 
final images on total integration time and atmospheric seeing are 
shown in Fig.|6]3, b. However, there are considerable deviations 
for individual tiles. F16Ks, for instance, turns out to be much 
shallower than expected based on its substantial total integra- 
tion time which is a consequence of the joined effect of worse- 
than-average seeing and the lower instrumental sensitivity in this 
period (see also Fig.|3j. The median of the depth of the final im- 
ages is 25.2 for J, and 24.7 both for H and A's, which is in good 
agreement with the original survey goals (cf. Sect. |2]l. Scaling 
the limiting magnitude to correct for the primary determining 
factors, namely integration time (fintg), FWHM, and the effec- 
tive variation of the instrumental zero point for this image, AZP, 
according to 

mii„ - 1 .25 log fintg + 2.5 log FWHM - AZP, (3) 



leaves a variation about the mean of between 0. 1 1 and 0. 14 mag 
RMS (Fig. [6];). This is caused by sky brightness fluctuations. 
There is the simple effect that the sky noise is higher when the 
sky is brighter - i. e., shot noise variations due to the fact that the 
sky is brighter on some nights than on others. Moreover, short- 
term fluctuations in the atmospheric sky brightness on a time 
scale similar to the integration time of an OB ultimately limit 
the accuracy with which the sky background can be subtracted 
during image data reduction. 

To quantify the depth of the survey as a whole, the effective 
area, that is the cumulative area distribution as a function of lim- 
iting magnitude, has been computed (Fig.|7ll. As the overlap of 
tiles is taken into account, the effective area corresponds with the 
final survey mosaics. One reads off that within 90% of the nom- 
inal survey area the limiting magnitude is 24.8, 24.3, and 24.1, 
within 75% it is 25.0, 24.5, and 24.4, and within 50% it is 25.1, 
24.7, and 24.7 for J, H, and Ks, respectively. 

4.3. Vaiidation of tlie pliotometric caiibration 

We have conducted several checks to verify the consistency of 
the photometric calibration of the final survey images starting 
with an internal check based on multiple detections, followed by 
the comparison with an independently reduced and calibrated 
subset of the same raw data, and, finally, by comparison with 
external NIR photometric data having been published in the lit- 
erature. 

At first, we have analyzed sources having multiple detections 
to verify the internal consistency of the photometric calibration 
of the final survey images. To this end, photometric differences 
of multiple detections of non-blended sources were inspected 
as shown in Fig. [Sh-ii. Total magnitudes have been used in or- 
der to minimize any possible effect due to seeing variations be- 
tween different tiles. We do not see any systematic discrepancies 
in this test, rather, at large, we found that the scatter appears 
to be compatible with the formal error bars. For J and H, the 
sources are located in the overlapping region of adjacent survey 
tiles, whereas for band, we took advantage of the "supple- 
mentary" survey fields F25 and F26 which mark a part of the 
survey area that is covered twice. Consequently, in case of J 
and H, the amount of scatter about the zero line is enhanced be- 
cause the sources lie towards the border of the image tiles where 
the exposure time of the jitter observations is effectively reduced 
with respect to the central part of the image, and, in addition, 
where possible flatfielding errors may be expected to increase. 
Regarding the two fields in K^, the larger scatter of F25 over F26 
is due to the fact that F25Ks is about 0.4 mag shallower than 
F26 (cf. Sect. 14. 2t . Generally, the overall scatter of measured 
fluxes with respect to true fluxes is expected to be smaller than 
the scatter shown here. 

Next, we have investigated whether the adopted strategy of 
data reduction and calibration procedures affect the final photo- 
metric results . To this end, we hav e used the J and photomet- 
ric data from ISaracco et al.l (1200 L) which is based on the same 
set of raw data that went into the production of tile F16 of the 
GOODS/ISAAC survey and allows to carry out a comparison 
that is unaffected by varying observing conditions or different 
instrumental characteristics. Here, total magnitudes were used 
in order to be compatible with the co mbination of correct ed 2" 
apertures and BEST apertures used bv lSaracco et al.l (1200 ll) . The 
results, displayed in Fig. [8^-/, show that both data sets are pho- 
tometrically compatible, except for very few outliers. 

Furthermore, we have checked our survey's photometry 
against external NIR catalogs of the CDF-S region which 
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have been obtained and calibrated independently from the 
GOODS/ISAAC data (Fig. |9]l. To this end, common sources 
were identified using a maximum pairing distance of 1" from 
which only isolated i. e. non-blended sources were selected. For 
the comparison with the Two Micron All Sky Survey (2MASS, 
ISkrutskie et al.ll2006l) . we also made use of the SOFI data set 
which had been used before to establish the survey's global 
photometric solution in order to extend the magnitude range 
for the comparison towards the bright end. In this case, just 
point sources were considered. Fluxes were measured in aper- 
tures of varying size (while maximizing the signal-to-noise ra- 
tio) and were corrected for aperture losses to a common aper- 
ture diameter of 10" which is virtually equivalent to a total 
magnitude given the typical PSF We applied the transforma- 
tion betwee n the 2MASS and LCO photometric system as de- 
termined by ICarpenteJ (l200lh . and used J - Jimass - 0.022, 
H - //2MASS = 0.012, and A's - A's 2mass = 0.015, whereas ne- 
glecting any color terms. The trend of measured source fluxes 
being apparently fainter than 2MASS fluxes at magnitudes faint- 
ward of 17.5 AB is the result of a Malmquist-like selection effect 
and is also present in the comparison with other data sets that are 
significantly shallower than GOODS/ISAAC. Taking this effect 
into account, the data exhibits no evidence for a bias in the pho- 
tometric calibration with respect to 2MASS. The scatter of the 
data is basically consistent with the formal flux errors which are 
clearly dominated by 2MASS. It is worth noting that consistent 
aperture corrections turned out to be crucial in this test. In fact, 
a lack thereof results in a significant bias at the ^0.1 mag level 
which is visible despite the significant scatter 

As we did not apply corrections for nonlinearity in the course 
of data reduction, neither for ISAAC nor for SOFI data, differ- 



ential nonlinearity may in principle affect any photometric com- 
parison. To address this issue we followed the same procedure 
as described in Sect l3.5.2l and discarded relatively bright, poten- 
tially nonlinear objects from the photometric analysis (indepen- 
dently in the ISAAC as well as in the SOFI data), thereby reduc- 
ing nonlinearity effects in our data to nominally ^ 0.01 mag and, 
thus, making them basically negligible with respect to the other 
sources of error. In terms of 2MASS the reference sources used 
here are relatively faint, which means that the photometric un- 
certainties are dominated by background noise and nonlinearity 
is expected to be insignificant. 

For the following comparisons with catalogs which are com- 
paratively deeper than 2MASS, we used SExtractor's auto- 
scaling magnitudes for formal consistency with the type of mag- 
nitudes published. However, photometric apertures were not 
matched for lack of detailed information about the apertures and 
the associated flux corrections for the reference data sets. Thus, a 
systematic effect depending on the data set is to be expected. The 
comparison with the EIS-DEEP data set in JHK^ is based on the 
most recent processing done by Olsen et aTl (l2006l) . The differ- 
ences of the AB corrections adopted in their and in our work are 
corrected for. It turns out that magnitudes from GOODS/ISAAC 
and EIS-DEEP are consistent modulo a constant offset of » 
0.1 mag for J and H, and certainly <0.1 mag for K^. The ap- 
parent scatter is larger than what the individual flux error esti- 
mates suggest - a finding that applies to this and the following 
comparisons as well. The magnitude diffe rences with respe ct to 
the //-band observations of the CDF-S by Mov et all (12003) are 
noticeably more dispersed while not exhibiting any systematic 
offset. T he comparison wi th the Las Campanas Infrared survey 
(LCIRS. IChen et al.ll2002h features a sHght trend with magni- 
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Fig. 9. Photometric comparison with respect to public NIR catalogs of the CDF-S region. Each pa nel corresponds to on e ref- 
erence data set in one filte r as indicated i n the label, namely, (a-c) Two Micron All S ky Survey (ISkrutskie et alJ 120061) . (d- 
f) ESO imaging survey (Olsen et alj 1200 6*). (g) //-ban d observations of t he CDF-S bv iMov et al] ^2003), (h) Las Campanas 
Infrared survey (Chenetal. 20021), and, (i) K20 survey dCknatti et al.l 12002') . The differences of magnitudes of sources as mea- 
sured on the GOODS/ISAAC images and as given in the respective reference catalogs are displayed as a function of the measured 
GOODS/ISAAC magnitude in the AB system. In case of 2MASS, dot symbols correspond to measurements performed directly on 
the calibrated ISAAC survey fields while square symbols denote measurements on the "auxiliary" SOFI images which were used to 
homogenize the photometric solution (see text). For the comparison with the other -much deeper- catalogs crosses and circles cor- 
respond to point-like and extended sources, respectively, with 1-sigma error estimates on the magnitude difference being indicated 
for point-sources if exceeding the size of the plot symbol. 



tude in such a way that bright sources (// ^ 18) appear brighter 
in GOODS/ISAAC while faint sources appear to be unbiased or 
even slightly fainte r than in the LCIRS . The magnitudes of the 
K20 galaxy survey dCimatti et al.ll2002 h appear to be brighter by 
a; 0. 1 mag than the respective GOODS/ISAAC magnitude. 

As a result, we do not find that photometric differences with 
respect to various independent data sets exhibit any coherent 
trend being attributable to the GOODS/ISAAC survey. Rather, 
it seems to be plausible to suppose that significant photomet- 
ric differences are mostly due to the photometric calibration, the 
exact definition of photometric apertures, the effect of seeing 
variations and, possibly, the lack of correction thereof of indi- 



vidual reference data sets. A mismatch of photometric apertures 
can make up a systematic difference of ~ 0.1 mag. For example, 
correcting the GOODS/ISAAC magnitudes for aperture losses 
by using total magnitudes completely eliminates the systematic 
offset with respect to K20. After all, the excellent photometric 
agreement with 2MASS in J and provides strong evidence 
that the global photometric error of the GOODS/ISAAC survey 
tiles is significantly below the 0.1 mag level. 
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Fig. 10. Performance of photometric estimators as inferred from image simulations. Ratio of recovered flux and input flux (magni- 
tude scale) for (a) Kron-like apertures (MAG AUTO), and (b) total, i. e. aperture-corrected, flux (MAG TOTAL). The binned average 
is plotted for J, H, and Ks with square, diamond, and circle symbols, respectively. The dotted line indicates the RMS scatter per bin 
for A^s- Panel (c) contrasts the actual flux errors computed in the RMS sense per bin (square, diamond, and circle symbols) with the 
flux error estimate according to Eq.|4](dashed, dash-dotted, and continuous line). The standard SExtractor flux error estimate for 
is shown as dotted line. 



5. Source catalogs 

5. 1 . Source detection and photometry 

We use SExtractor dBertin & ArnoutsI 1 996|) for source detection 
adopting a spatial filtering with a 2-dimensional Gaussian of 2.5 
pixels FWHM to match with the PSF of the best seeing images. 
A detection threshold of 1.35 times the local background RMS 
has been applied with a minimum detection area of one pixel. 
This choice corresponds to a 5cr significance once the spatial 
filtering is taken into account. The choice of detection thresh- 
old and minimum detection area has been gauged based on the 
analysis of catalog completeness and contamination for a grid of 
detection parameters (cf. Sect. 16.1b . For the given spatial filter- 
ing, we have tested in particular if a minimum detection area of 
more than one pixel is favorable. However, it turned out that this 
does not increase the catalog completeness at a given contami- 
nation level, neither for point-sources nor for extended sources. 
The background map has been estimated in square cells of 200 
pixels on a side using a 5 x 5 cells median filtering. For photom- 
etry a local background annulus in square shape with 20 pixels 
width was used. 

We have adopted the prescription of iLabbe et"an (I2003L 
Sect. 5.2) to define the aperture-corrected total flux and a flux 
from which colors are computed based on the Ks band aper- 
tures. For isolated sources the total flux measurement is based 
on SExtract or's au to-scaling elliptical apertures (MAG AUTO) 
inspired by iKronI (|1980) corrected for the flux loss due to the 
finite aperture size. For the parameters which control the size 
of the auto-scaling apertures, we have used the default values, 
that is A: = 2.5 for the scaling parameter and R^in = 3.5 for the 
lower bound of the aperture size (in isophotal units). To apply 
the image-specific aperture corrections as given in Sect. 14.11 we 
have converted the elliptical aperture area into the correspond- 
ing circular aperture's diameter In order to minimize the mutual 
flux contamination for blended sources, in this case the total flux 
measurement is based on a reduced circular aperture whose area 
is half the isophotal area. To prevent too small and too large aper- 



tures that are more error-prone, the aperture diameter is bound 
by 0.7" and 2.0", respectively. The "color flux" is isophotal (as 
defined in K^) if the source is isolated and if the isophotal area 
ranges between the size of a circular aperture of 0.7" and 2.0" 
diameter For sources whose isophotal areas exceed these limits, 
the aperture flux at 0.7" or 2.0" is assigned instead. For blended 
sources, the color flux is equally defined as for the total flux, 
based on the flux measured within a "reduced" circular aperture. 

We computed flux errors taking into account the background 
fluctuations, both for the source aperture and for the background 
region, and the contribution due to the discrete nature of the 
detector electrons. To this end, the noise model established in 
Sect. 14. 21 was evaluated at the spatial scales corresponding to the 
aperture size. A, and the size of the background region, B, re- 
spectively. With the discreteness noise which has been estimated 
based on the local gain, G, the flux error is given by 

A/= (4) 

We investigated the performance of flux and flux error mea- 
surements by means of simulated images, that is actual survey 
tiles into which stellar sources of known flux had been embed- 
ded (cf. Sect. |6X2] |. Fig. [TOj? shows that SExtractor's MAG 
AUTO measurement suffers a systematic flux loss of ^0.1 mag 
at 20 mag that increases towards fainter fluxes, amounting s; 
0.2 mag (at 24 mag) for isolated sources throughout all bands. 
The apparent reduction of the average flux loss faintward of 
24.5 mag is due to the selection effect which renders the distri- 
bution to be significantly skewed towards larger measured flux 
(Malmquist bias), and since this bias is determined by the signal- 
to-noise ratio, the effect for the J band is less pronounced than 
for H and K^. Fig.fTOb demonstrates that the measured total flux, 
denoted MAG TOTAL, is unbiased. 

Fig, [rob illustrates that the flux error computed according to 
Eq.|4]is able to represent the true flux errors quite accurately over 
a wide range of fluxes. In contrast, SExtractor's standard flux 
error (MAGERR AUTO) systematically underestimates the true 
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Fig. 11. NIR color-magnitude diagram (main plot) and J - Ks 
color-redshift distribution (inset). The color-coding is accord- 
ing to spectroscopic redshift using seven intervals between 0, 
0.5, 0.67, 0.82, 1.04, 1.23, 2, and 6, as shown in the inset. 
Sources without spectroscopic redshift are shown as black dots. 
Sources which appear point-like in the HST/ACS z-band image 
are marked by crosses. 



flux errors by at least a factor of ~ 1 .5 (in mag) for all fluxes. This 
discrepancy underlines that it is essential to take into account 
the actual local background fluctuations for realistic flux error 
computations instead of simply scaling the pixel-to-pixel vari- 
ance. In detail, the estimated flux error is in fairly good agree- 
ment faintward of a; 22 mag. The apparent overestimation of the 
flux errors for the faintest fluxes (S 24 mag, especially for H and 
Ks) is probably an artefact due to the strong selection effect in 
this regime. Towards bright fluxes, errors are increasingly un- 
derestimated, for instance, by about a factor of two at ~ 20, 
because systematic effects that are beyond the simple statistical 
error model, such as source blending, start to become dominant. 



5.2. The Color Catalog 

To create the selected NIR color catalog, we ran SExtractor 
for each of the 78 tiles in double image mode with the band 
mosaic for detection and with the respective tile for measure- 
ment using the photometric scheme described in the previous 
section. This approach guarantees consistent and accurate flux 
measurements by using unique apertures defined in /Tj for all fil- 
ters and by allowing to take into account the variation of the ef- 
fective PSF and the sky background noise from image to image. 
The detections from all tiles were merged to obtain the final cat- 
alog of unique sources. In case of sources which were detected 
multiple times in the same band - in regions of tile overlap - the 
best signal-to-noise measurement evaluated based on the local 
limiting magnitude was kept (cf. Sect. 14.2b . Eventually, differ- 
ential aperture corrections were applied to the measured colors 
using Table |3] The resulting catalog reports J - and H - 
and the corresponding Icr uncertainties for 7079 sources which 
have solid, i. e. >5cr, measurements in all three bands. Note that 
corrections for galactic extinction were not applied. The cata- 



log data can be downloaded directly from the ESO archive by 
navigating to the URL given in Sect. [T] 

The catalog columns are as follows. Columns no. 2 to 1 1 are 
directly adopted from SExtractor run on the images. 

(I) NUMBER - sequential source number for reference. 

(2, 3) ALPHAJ2000, DELTAJ2000 - Source position R.A., 
Declination (J2000). 

(4) ISOAREAFJMAGE - isophotal area in pixel units. 

(5) KRON_RADIUS - Scaling of the automatic elliptical aper- 
ture inspired by Kron ( 1980). 

(6) FWHM JMAGE - FWHM measured in pixels units. 

(7) ELLIPTICITY - Source ellipticity. 

(8-10) A.WORLD, B.WORLD, THETA.WORLD - ElHptical 
shape parameters: semi-major axis, semi-minor axis and ori- 
entation (in degree). 

(II) FLAGS - SExtractor flags, possibly combined in a logical 
OR fashion. 1 - object has close neighbors, 2 - object was 
blended, 4 - object is saturated, 16 - aperture incomplete. 

(12, 13) J_MAG.COLOR, J.MAGERR_COLOR - J magni- 
tude and associated error from which colors are computed. 

(14-17) As before but for H and K^. 

(18, 19) K.MAG_TOTAL, K.MAGERR.TOTAL - Total 
magnitude assuming a point source profile and associated er- 
ror 

(20, 21) J_K, J_K_ERR - J - color and associated error. 
Differentially corrected for aperture losses. 

(22, 23) H_K, H_K_ERR - H - K, color and associated eiTor. 

(24) J_FRAME - Original survey tile from which the J mea- 
surement was extracted. This allows to trace back to the orig- 
inal data if needed, for instance for the creation of image 
cut-outs. 

(25, 26) As before but for H and K^. 

Given that aperture corrections derived for point sources are ap- 
plied to mostly faint galaxies, in general, this is likely to yield an 
underestimate of the total flux (column 18), since the light pro- 
file for galaxies will be more extended than that of point sources. 
The amount of this effect will be quantified as a function of the 
extent of galaxies in Sect. |6.1.3] using image simulations. 

5.3. Color distributions 

To validate the observational data the color distribution is exam- 
ined and put in perspective with redshift data taken from the lit- 
erature. To this end spectroscopic redshif ts were compil e d from 
a number of published data sets, namely 'Cristiani et al.' (i2000l). 
Croometal. (2001), Le Fevre et al. (2004), Strolg er et al.1 
(2004), Szokolvetal. (2004), van der Wei et al.l " (|2004. 
Mignoli et al. (2005), Vanzella et al.. (2008), and Popesso et al.l 

Fig. [TT| shows the J - color-magnitude of /Tj -selected 
sources in the CDF-S. Owing to their bluer J - color, the 
locus of stellar sources appears sharply offset from the bulk of 
redshifted galaxies which, in contrast, populate a large range of 
colors. In sufficiently deep observations red NIR color selection, 
J - Ks > 1.34 (corresponding to - K^ ^ 2.3 in the Vega 
system), is an efficien t method to find high-redshift, i. e. z ^ 2, 
galaxies (iFranx et al.ll2003b . In our data, very red objects having 
J - Ks > 1.34 are found in the Ks range between 21 and 24.5 
peaking at Ks ~ 22.1 . However, we note that to date there are 
very few spectroscopic redshifts for distant red galaxies in this 
field (see, e. g., Kriek et al. 2006). 

The NIR color-color diagram is shown in Fig. [12] with color 
tracks of four Hubble-types (El, Sbc, Scd, Im) and two starburst 
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Fig. 12. Observed NIR color distribution overlaid with tracks of El, Sbc, Scd (left panel), and Im, SB3, and SB2 (right panel) 
galaxy spectral templates for a redshift range between 0.01 and 3 (as indicated along the curves). The color-coding according to 
spectroscopic redshift has been adopted from Fig.[TT] Average color errors ( 1 cr) being representative for sources with = 24, 23, 22 
are displayed top left. Sources which appear point-Uke in the HST/ACS z-band image are marked by crosses. 



ga laxies (SB2, SB3) ov erlaid for which the spectral templates 
of 'Coleman et al.' (1980) and two templates of starburst galaxies 
from Kinnev et al. ( 1996) were redshifted to between 0.1 and 3 
and folded with the actual instrumental throughput. Part of the 
BPZ software package of B enftez (2000) has been applied for 
this calculation. The smooth tracks of the El, Sbc, and Scd-types, 
originate from the sampling of their uniform spectral energy dis- 
tribution which is dominated by the absorption-modulated con- 
tinuum. In contrast, the young stellar populations of the starburst 
types (SB3, SB2), and to some extend also the Im-type, with 
their strong spectral emission features result in more wrinkled 
tracks. 

The noticeable loop between redshift 0.3 and 0.5, at which 
the H ~ Ks color gets bluer with increasing redshift is typical 
for an old stellar population and therefore most prominent for 
the El, Sbc and Scd types. It occurs when the restframe band 
samples the strong absorption bands shortward of 1 .5 fim mainly 
due to TiO. For an old stellar population seen at redshifts be- 
yond one, J - Ks is strongly increasing because the restframe J 
flux is continuously fading due to the strong stellar atmospheric 
absorptions mainly induced by metals like Mgl, Fel, Cal etc. 
Towards higher redshift above ~ 1.7 this trend becomes even 
more exaggerated as the 4000- A/Balmer break starts to enter the 
J band. The redshift dependence of the NIR colors for redshifts 
above 1 .5 illustrates the merit of NIR data for the determination 
of accurate photometric redshifts when the characteristic absorp- 
tion features have moved out of the optical bands. A very young 
stellar population seen at high redshift can not reproduce the 
very red NIR color as seen in the observations. Therefore, the 
detected very red objects can be interpreted as distant galaxies 
hosting an evolved stellar population. Otherwise, heavily dust- 



absorbed star burst galaxies could give rise to such very red col- 
ors. 

6. Analysis 

6.1. Completeness and contamination 
6.1.1. Empirical estimation 

Under the ESO programme 73.A-0764, lead by I. Labbe, a very 
deep 21 -hours ISAAC integration of the Hubble Ultra Deep 
Field (UDF) inside of the CDF-S region was obtained. The 
pointing is centered on F14 whereas the instrument was rotated 
by ca. 45 degrees to optimally match the Hubble UDF observa- 
tions (Fig.lTJ. We have reduced the raw data, which are publicly 
available from the ESO science archive, using similar procedures 
as for the whole GOODS/ISAAC survey data reduction (Sect.[3j. 
Furthermore, in order to improve the background subtraction for 
faint sources, we have applied a third pass in the OB process- 
ing using a source mask generated from the co-added OBs. The 
photometric zero points were bootstrapped from the FlAKf^ im- 
age. We have co-added the calibrated OBs including the FIAK^ 
GOODS/ISAAC data resulting in the very deep image coving 
5.40 arcmin^ of the UDF region with a total integration time of 
91 182 sec, an eff'ective PSF of 0.36" FWHM, and a 5cr depth for 
point sources of 25.60 mag which is 0.67 mag deeper than the 
respective survey image (see also Fig.jSJ?). 

We have utilized this image as a reference to empiri- 
cally study the completeness and contamination, i. e. the frac- 
tional number of missed sources and spurious detections, in the 
GOODS tile F14 in the band. As this is a typical tile in terms 
of total exposure time, image quality etc. (see Sect. 14. It , the re- 
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Fig. 13. Central 1' X 1' cut-out of (a) the 6-hours GOODS/ISAAC survey field F14 in K„ and (b) the 25.3-hours co-added 
UDF/ISAAC image, both displayed at linear scaling. The effective PSF has 0.46" and 0.36" FWHM, and the 5cr depth for 
point sources is 24.93 and 25.60 mag, respectively. 



suits can be considered to be representative for the whole sur- 
vey in Ks band, and - to a minor degree - for the J and H 
bands. To identify spurious detections and missed sources, we 
have cross-correlated the source catalogs of the two images us- 
ing a 1" matching distance. Fig. [14] displays the source counts 
aggregated within 0.35-mag bins according to the total flux and 
the resulting completeness and contamination levels with respect 
to the reference image. The relative completeness was calculated 
from the ratio of the number of cross-identified sources and the 
total number of sources in the reference image both counted with 
respect to the total flux measured in the reference image. Hence, 
the resulting numbers purely represent the detection efficiency 
without being affected by either spurious detections or the mi- 
gration effect due to flux measurement errors (Eddington bias). 
Counts faintward of - 24.8 start to become significantly in- 
complete, dropping to a completeness of 46% at = 25.15. 
At this flux level, the relative contamination which is defined 
as the ratio of counts of spurious sources over counts of true 
sources, both based on the total flux measured on F14, amounts 
to 26%. Despite the increasing statistical uncertainty towards 
fainter fluxes, it is evident that the contamination becomes a 
dominant effect beyond the 50% completeness magnitude. 

6.1.2. Image simulations 

We have executed extensive image simulations in order to char- 
acterize the source detection efficiency i. e. the completeness 
level of the resulting source catalogs as a function of source flux 
for each tile. 

The simulations were set up as follows: for each survey tile 
image, first, an isolated unsaturated star was selected and its 
7.5" X 7.5" image was cut out to create a representative template 
that is suitable to quantify the completeness for point sources in 
that tile. In order to estimate the completeness for faint galax- 
ies, we have smoothed the stellar images using a Gaussian ker- 



nel of 0.35" FWHM, which has been adjusted so that the dis- 
tribution of central surface brightness of simulated faint objects 
(~24 mag AB) matches the observational data. Although this is 
certainly not sufficient to simulate the observed range of galaxy 
types, for example low surface brightness galaxies, this approach 
allows the completeness to be estimated at first order without 
resorting to extra morphological parameters. A more accurate 
treatment of the statistical incompleteness taking into account a 
proper mix of morphological types is beyond the scope of this 
work but will be addressed in a forthcoming paper which will 
present a comparison with physical models of galaxy evolution. 
Then, the template image was dimmed to a flux in the magni- 
tude range between 20 and 27 which was drawn at random from 
a power-law distribution whose slope of 0.26 (per magnitude) 
resembles the typical number counts in the field. That followed, 
the dimmed simulation template was co-added to the image at 
a random position. Per image, 500 realizations were simulated, 
each of which with 100 simulated sources. This corresponds to 
a relatively small fraction of artificial sources to real sources in 
the images of 1:8.5 on average so that the extra crowding and 
the statistical incompleteness due to blending induced thereby is 
negligible. Source detection and photometry was carried out on 
the simulated images in exactly the same way as for the obser- 
vational data using identical detection parameters. For each tile, 
we have delimited the analysis to the central 75% of the image 
area within which the exposure time is approximately uniform. 
Fig. [T4b shows that the completeness resulting from the sim- 
ulations is consistent with the independent empirical estimate. 
This finding demonstrates the viability of the adopted simula- 
tion approach in spite of its inherent simplicity of how galaxies 
are modeled. 

Based on the image simulations, we have characterized each 
tile in terms of the completeness function for point sources 
and for galaxies. On average the galaxy completeness is about 
0.3 mag brighter than the point source completeness. The strong 
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Fig. 14. (a) Ks source counts in the GOODS field F14 with 
respect to the 0.67-mag deeper UDF/ISAAC data using ex- 
actly the same region. Plot symbols are as follows: all sources 
detected in the UDF/ISAAC image - squares, true positives, 
that is, objects cross-identified between F14 and UDF/ISAAC 
- dashed line, missed sources, that is, objects that are only de- 
tected in UDF/ISAAC - dotted line, and spurious detections, 
that is, objects that are detected in the GOODS F14 data but 
not in UDF/ISAAC - triangles. Error bars are Poissonian. The 
small arrows on the top indicate the 5cr limiting magnitude of 
point sources for the two images, (b) Resulting relative source 
completeness (squares) and relative contamination (triangles) in 
comparison to the completeness estimate based on image simu- 
lations (thick line, red). 



correlation of mso, defined as the total flux for which the point 
source completeness drops to 50%, with limiting magnitude can 
be parameterized by niso ~ mii,n(5cr) + 0.25. 

6.1.3. Surface brightness selection 

Surface brigh tness selection eff ects are inevitable in any galaxy 
sample (e. g. iLillv et aT]|1995h as well as the loss of flux for 
extended sources. To quantify the two effects, first, spheroids 
and disks were modelled using deVaucouleurs and exponential 
profiles, respectively, and were embedded into the image of 
field F14. As a surface brightness measure, we adopt the surface 
brightness as inferred from the detected peak flux of each source 
image. Fig.fTSb and c show that the surface brightness serves as 




Fig. 15. Completeness, peak surface brightness and flux loss for 
spheroids (top panels) and disks (bottom panels) as inferred from 
image simulations for the GOODS field F14. denotes the ef- 
fective radius for spheroids and h the scale length of disks. Grey 
shading displays the completeness in levels between 10% and 
90% in each plot. Peak surface brightness (in mag arcsec"^) is 
overlaid on the left hand side (a, c), and flux loss (in steps of 
0. 1 mag starting from 0. 1 mag) is shown on the right hand side 
(b, d). 



primary selection parameter, almost independently of total flux, 
source extent, and type of profile. In this case, the 90% com- 
pleteness limit corresponds to a peak surface brightness of 24.0 
mag arcsec"^, while the completeness drops below 50% for 24.5 
mag arcsec"^. For other images, however, the numeric values are 
expected to vary according to image depth. Fig.flSt' and d illus- 
trate the flux loss (for MAG TOTAL) featuring a strong increase 
with increasing source extent. Part of the substantial flux loss of 
large disks is due to source blending which inhibits the proper 
measurement of the total flux. 

The surface brightness selection effect at 24.0 mag arcsec"^ 
also becomes apparent in the observed distribution of galaxies 
(Fig.fTSIl. Thanks to the deep integration of the UDF region, the 
properties of the sources being missed in the GOODS survey can 
be directly examined. In the two faintest bins in which galaxy 
counts are being accumulated, marked by vertical dashed lines, 
centered at = 24.5 and 24.75, the bias against low surface 
brightness objects is evident. Consequently, the resulting total 
counts in this range may be biased, but it is difficult to quantify 
the net effect. For objects brighter than ~ 24.4 mag, however, 
there is apparently no significant selection effect. 

To further examine the selection of faint galaxies depend- 
ing on their physical and morphological properties, we extracted 
images of moderately distant objects (z between 0.15 and 0.6), 
scaled them according to redshift taking into account cosmologi- 
cal dimming and the K-correction, and computed surface bright- 
ness and flux as if they were actually observed. Fig. [T6] displays 
the resulting tracks of galaxies with the following characteris- 
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Fig. 16. Peak surface brightness as a function of K^; magnitude 
as defined in the text. Detections in GOODS field F14 (solid 
dots) opposed to missed sources (open circles) that were identi- 
fied in the deeper UDF/ISAAC image. The few stars in the field 
are marked by cross symbols (cyan color). The insert shows the 
equivalent diagram for the GOODS survey as a whole with point 
sources being highlighted (cyan). Overlaid tracks, labelled from 
a to d, correspond to individual galaxies as if they were observed 
at the redshifts denoted along the curves. 



tics: (a) 2.1 M* spheroidal, (b) 0.1 M* disk, exponentially domi- 
nated profile, (c) 0.07M*, low-surface brightness, sub-dominant 
central bulge otherwise irregular morphology, (d) 2.8M* nearly 
face-on spiral, disk-dominated with central bulge. The luminous 
spheroidal, a, which has the most centrally concentrated profile 
of the sample, is least affected by surface brightness selection. 
In contrast, the spiral galaxy, seen nearly face-on, d, appears ca. 
1 .3 mag fainter in surface brightness which would be hardly de- 
tectable beyond z = 2.5 although the luminosities of the two 
galaxies are similar. 

6.2. Number counts 

Differential galaxy number counts were independently deter- 
mined for each of the three NIR bands as follows. We have 
counted the number of sources in 0.25-mag bins with respect 
to the total flux (as defined in Sect. 15. It whereas stellar sources 
had been excluded on the basis of their point-like appearance in 
the GOODS/ACS z-band mosaic (see below). First, each tile was 
processed separately and the raw counts were corrected for in- 
completeness according to the completeness function for galax- 
ies which had been established using simulations (Sect. l6Tt . At 
this point the central part of each tile, 75% in terms of area, 
was used over which the exposure is approximately constant. 
Then, for each magnitude bin the contributions from individual 
tiles having at least 50% point source completeness were com- 
bined. In this way the variance of the effective area as a function 
of depth was accounted for. The correction factor is eff'ectively 
bound and, consequently, only the deepest tiles contribute to the 



faintest bins (which are J and Kg - 25.25 and H - 25.0). The 
actually applied completeness correction factors are about 1.6 
for the faintest GOODS/CDF-S data points in all three bands, 
and 1.7 for the data point at = 25.25 which was extracted 
from the UDF/ISAAC image. The resulting diff'erential counts 
are displayed in Fig. [17] and reported in Table |4]where cumula- 
tive counts, A^(< m), are quoted in addition. 

As identification criterion for stars, a cut in R^o was used 
in conjunction with the flux Umit zab < 24.5 (MAG_AUTO) 
to restrict the selection to the flux range where the stellarity 
criterion is able to clearly identify the stellar locus. Then, it 
was verified that all sources have J - Ks < 0, which is a very 
good criterion for the whole range of stellar types except for 
the very cool but c omparatively rare L and T dwarf types (e. g. 
iHawlev et al.ll2002b . leading to the removal of 5 objects most of 
which have been independently spectroscopically identified as 
QSO (cf. Fig. fTTT i. The number counts of the resulting sample 
of 295 stars are perfectly concor dant with the Galaxy model as 
simulated bv lGirardi et al.l (l2005h (Fig.[T7] right panel). Beyond 

~ 24, the sample becomes incomplete due to the imposed flux 
limit in z band, as expected. Furthermore, using the angular cor- 
relation function, we have checked that the magnitude-limited 
stellar sample (17 < < 23) is uncorrected over the full range 
of angular scales probed by the survey (0 ~ 0.0005 . . .0.1°), 
which, besides, confirms the photometric homogeneity of the 
data. Altogether, the tests provide no evidence for a significant 
contamination of the stellar sample from e. g. compact galaxies. 
Thus, the resulting galaxy number counts are simply aff'ected by 
the incomplete subtraction of stellar sources beyond ~ 24. 
The amplitude of this eff'ect, however, is just about 2% which is 
less than the size of the plot symbols in Fig.fTT] 

The faint-end galaxy number counts (m > 22) follow an ex- 
ponential with logarithmic slope of a - (0.262 + 0.007) mag"' 
for J and, marginally shallower, a - (0.254 ± 0.008) mag"', for 
H and K^. The best-fitting amplitudes (in a least squares sense) 
are 3.31, 5.77, and 6.40 x 10"- mag"' deg"^ for J, H, and K„ 
respectively. Down to the faintest fluxes we do not find an indi- 
cation for a significant change in the slope of the number counts. 
According to the size of the survey area which exceeds previ- 
ous surveys of similar depth by about an order of magnitude, 
the statistical errors turn out to be almost negligible over a large 
range of fluxes. At the faintest fluxes, the statistical errors of 
the data points must slightly increase because of the decrease of 
the effective area. Generally, the results in the low-flux regime 
may be systematically affected by uncertainties in the incom- 
pleteness correction and possible contamination due to spurious 
detections. To minimize the possible impact, we have limited the 
analysis to ^ 5cr significant sources. As a consequence of this 
choice, correction factors are quite moderate and the contamina- 
tion level remains low (~25%). 

Comparing our results with previously published data of 
deep surveys selected from the literature, we find a good agree- 
ment in general (Fig. [T7] l. However, closer examination reveals 
some discrepancies much in the same way as discrepancies 
have been found amongst p revious data se ts (e. g. ISar acco et aTl 
120011: iMetcalfe et aLl l2006'). For instance, 'Sara cco et all (|200"TI) 
reported, based on the deep imaging of the F16 field in the CDF- 
S region, comparatively steep 7-band number counts with sys- 
tematically more faint {J > 24) galaxies than we have found. 
By contr ast the number coun ts resulting from the Subaru Deep 
Survey ( Maihara et al.ll200TI) are consistent with our results in 
the 7-band (and similarly in the A's-band), likewise the raw 
number counts in the Hubble Deep Fi eld-North obtained wit h 
HST/NICMOS using the FllOW filter ([Thompson et al.|[T999h . 
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Fig. 17. Differential galaxy number counts in the GOODS/CDF-S region in the J, H, and Kg bands (left, middle, and right panel, 
respectively) corrected for incompleteness and displayed with Poissonian ( 1 cr) error bars. The two faintest data points in (asterisk 
symbols) were inferred from the very deep UDF/ISAAC image. The dashed lines indicate the faint-end exponential fits with log- 
arithmic slopes between 0.24 and 0.27 mag"' (see text). Results from other NIR surveys selected from the literature are displayed 
for comparison. Generally, incompleteness-corrected data is shown except for th e HST/NICMOS resu lts for which incompleteness 
at the depth of the GOODS/ ISAAC data is negligib le. Note that the data labeled lSaracco et al. l (l200TI) refers to the results obtained 
from the CDF-S, the data of Thomps on et al] J1999I) shown here is limited to 7 > 24 for which the relatively small field provides a 
sufficient absolute number of counts. The apparent bump around Ks ~ 20 in the data of Forster Schreiber et al. (2006) is owing to 
the contribution from the MS 1054-03 clus ter of galaxies on which the su rvey field is centered and shoul d not be seen as disc repancy 
with the other data sets. The data from lForster Schreiber et all (12006 1) HDF-S refers to the results of iLabbe et al.l (l200 3l) but it is 
apparently not identical, however, as to the comparison with our results this incongruity see ms to be irrelevant. In the right panel, the 
stellar A's-counts in the GOODS/CDF-S (circles) are shown in contrast to the simulation o f iGirardi et all (l2005l) based on a model 
of the Galaxy (dotted line). 



In the //-band the results of two imaging surveys which have 
been c ond ucted in the CDF-S region previously, by iMov et al] 
(l2003h and ICimatti et a l. (2002), are fully compatible with our 
findings. Also the number count s resulting from the completely 
independent study by Metcalfe et al.l ( 12006) based on imaging 
of the 7' X 7' William Herschel Deep Field using the Omega 
Prime camera on the 3.5-m Calar Alto telescope, Spain, are 
perfectly in agreement with our data. The extrem ely deep H- 
band number counts reported by iMetcalfe et al.l (|2006) based 
on observations of the Hubble Deep Field-South (HDF-S) us- 
ing NICMOS seem to be slightly in excess with respect to our 
findings, however, taking into account the small size of the HDF- 
S, 0.95 X 0.95 arcmin^, the two results are consistent within 
the Icr Poissonian errors. Our /Tj-band results are compatible 
with the UKIDSS UDS, that is the ulti-a-deep survey of the 
suite of infrared de ep sky survey s using the United Kingdom 
Infrared Telescope (iHartlev et al . 2008), the faint infrared ex- 
tragalact ic survey (FIRES) in the MS 1054-03 field and i n the 
HDF-S (iForster Schreiber et al.1 l2006t iLabbe et all l2003h . the 



K20 survey [ Cimatti et al. I (12002), the TiTs-band observations of 
'Saracco et al. (2001) in the CDF-S and the Subaru Deep Survey 
(Maihara et al. 2001) once the formal uncertainties are taken into 
account. Concerning the faint end of the counts, though indi- 
vidual data points are formally consistent within the error bars, 
the slope that we find is steeper than w hat has been r eport ed for 
the FIRES survey, namely a a; 0.15 bv'La bbe etaTI (12003 ). and 
a = 0.20 by Forster Schreiber et al. (2006). However, the com- 
parison of number counts results must be exercised with caution. 
At faint flux levels, close to the limiting magnitude of the sur- 
veys, the statistical corrections and their difficult to assess un- 
certainties usually dominate the error budget. In general, any 
photometric mismatch, e. g. due to discrepancies in the pass- 
band, the details of the flux measurement scheme, or the correc- 
tion to total flux, may result in differences of the number counts. 
Furthermore, the presence of cosmic large-scale structure (LSS) 
causes a variation of number counts from field to field in excess 
to the usually quoted fluctuation according to Poissonian statis- 
tics. 
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To quantify the effect of cosmic variance on our results, we 
have assumed an angular correlation of power-law form with the 
canonical index, i. e. w{0) = Afl""^, and a correlatio n amplitude 
based on the measured angular clustering (see, e. g.. lDaddi et alj 
12003 1). Taking into account the actual survey area, we find the 
fractional cosmic variance of the cumulative counts for m < 22, 
23, and 24 mag (AB) to be 3.4, 3.2, and 3.0% (Icr), respec- 
tively, in addition to the respective Poissonian fluctuations of 
2.5, 1.8, and 1.7%. Put differently, cosmic variance increases 
the Poissonian errors on the (cumulative) number counts by fac- 
tors between ~ 1.7 and 2.1. According to the power-law model, 
Poissonian fluctuations decrease more rapidly with increasing 
survey area (Q.) than the LSS fluctuations, namely oc Q^"^ vs. 
oc Q Therefore, the contribution of cosmic variance to the to- 
tal error budget will be significant also for next generation deep 
surveys. For number counts down to 24 mag in a 1° x 1° field, 
for instance, one expects about 0.34% Poissonian and 1.6% LSS 
fluctuations. 



7. Summary 

We presented the deep VLT/ISAAC NIR imaging survey in the 
GOODS-South field. The survey was conducted as part of the 
ESO/GOODS program and the resulting data products were re- 
leased to the public via the ESO science archive. The survey 
covers an area of 173 arcmin^ in J and K^, and 160 arcmin^ in 
the H band. The "median" survey depth is 25.1 in J and 24.7 
both in H and A's in terms of the 5cr total limiting AB magnitude 
for point sources. The excellent image quality is characterized 
by a PSF with FWHM between 0.34" and 0.65". The asti'omet- 
ric calibration is accurate to ~ 0.06" RMS over the whole field. 
The absolute astrometric accuracy is limited by the accuracy 
of the GSC 2 on which the astrometry is ultimately based. The 
photometric calibration was verified in various ways including 
tests against independent NIR photometric measurements such 
as 2MASS. The internal accuracy of the photometric calibration 
is < 0.017 mag (Icr). The overall photometric accuracy, includ- 
ing systematic effects, is up to ~ 0.05 mag for all three bands. 
We described the survey layout, observations, data reduction and 
calibration, and the quality control process. The data products 
were characterized in terms of PSF width, curve of growth, and 
limiting depth. A /Ts-selected catalog of 7079 objects sources 
with JHKs photometry was defined for which color distributions 
were presented with respect to available redshift data. 

We characterized the properties of extracted source catalogs 
in terms of their fractional completeness and contamination as 
a function of magnitude. For this purpose additional 21 hours 
of Ks imaging data of the Hubble Ultra Deep field, which are 
available from the ESO archive, were processed and the result- 
ing image, which is ~0.7 mag deeper than the respective survey 
image, was analyzed. These results were combined with exten- 
sive image simulations. For the number counts of faint galax- 
ies we found very similar logarithmic slopes in all three bands, 
ranging at the faint end between 0.24 and 0.27 mag"' without an 
indication for a change in the slope. Our new results are largely 
compatible with the ones from previous surveys. However, the 
GOODS/ISAAC data establish the surface density of faint galax- 
ies in the universe with unprecedented accuracy owing to its 
unique combination of survey area and depth which is unrivaled 
by any other NIR survey. In fact, this survey represents the deep- 
est NIR imaging data set to date covering a contiguous region 
well above 100 arcmin^. In the H band, there is to our knowledge 
no ground based survey paralleling GOODS/ISAAC in depth. 
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Table 2. Properties of the final survey tiles. 
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Table 2. continued. 
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Table 3. Aperture corrections of the final stirvey tiles for aperture diameters be- 
tween 0.7"and 8.38"assuming point source profiles. 
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0.116 
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0.143 
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0.128 
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0.115 
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0.124 
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0.131 
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0.096 
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0.093 
0.057 
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0.102 
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0.081 
0.078 
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0.080 
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0.072 
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0.121 
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0.117 
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0.092 

0.071 
0.077 
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0.069 
0.085 
0.062 
0.072 
0.077 
0.050 
0.069 
0.052 
0.065 
0.064 
0.049 
0.036 
0.072 
0.041 
0.051 
0.070 
0.070 
0.064 
0.072 
0.053 
0.052 
0.065 
0.059 
0.057 
0.066 
0.059 
0.049 
0.069 
0.053 
0.042 
0.069 
0.089 
0.065 
0.056 
0.065 
0.083 
0.076 
0.085 



0.067 
0.051 
0.055 
0.044 
0.051 
0.063 
0.043 
0.054 
0.058 
0.035 
0.050 
0.038 
0.045 
0.046 
0.034 
0.025 
0.054 
0.030 
0.036 
0.053 
0.052 
0.044 
0.052 
0.036 
0.037 
0.047 
0.043 
0.042 
0.049 
0.040 
0.033 
0.052 
0.038 
0.029 
0.051 
0.064 
0.049 
0.043 
0.046 
0.062 
0.055 
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0.047 

0.036 
0.038 
0.031 
0.037 
0.046 
0.031 
0.039 
0.042 
0.024 
0.037 
0.025 
0.031 
0.032 
0.022 
0.017 
0.036 
0.021 
0.026 
0.038 
0.037 
0.032 
0.034 
0.025 
0.026 
0.034 
0.030 
0.029 
0.035 
0.028 
0.024 
0.037 
0.026 
0.019 
0.037 
0.045 
0.034 
0.031 
0.032 
0.047 
0.038 
0.043 



0.031 

0.024 
0.027 
0.018 
0.025 
0.031 
0.020 
0.027 
0.029 
0.016 
0.024 
0.017 
0.020 
0.021 
0.015 
0.011 
0.024 
0.014 
0.017 
0.026 
0.026 
0.022 
0.023 
0.016 
0.017 
0.023 
0.019 
0.019 
0.024 
0.019 
0.016 
0.024 
0.017 
0.013 
0.024 
0.030 
0.023 
0.023 
0.021 
0.032 
0.026 
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0.019 

0.015 
0.016 
0.011 
0.015 
0.017 
0.013 
0.016 
0.018 
0.010 
0.016 
0.010 
0.013 
0.014 
0.009 
0.007 
0.015 
0.009 
0.011 
0.017 
0.017 
0.014 
0.014 
0.010 
0.011 
0.014 
0.012 
0.012 
0.016 
0.012 
0.010 
0.016 
0.011 
0.008 
0.015 
0.019 
0.014 
0.014 
0.013 
0.020 
0.017 
0.019 



0.010 

0.009 
0.009 
0.006 
0.008 
0.010 
0.007 
0.009 
0.010 
0.005 
0.009 
0.006 
0.007 
0.008 
0.005 
0.004 
0.009 
0.005 
0.006 
0.010 
0.010 
0.008 
0.008 
0.006 
0.006 
0.008 
0.007 
0.007 
0.009 
0.007 
0.006 
0.009 
0.006 
0.004 
0.009 
0.010 
0.008 
0.007 
0.007 
0.012 
0.010 
0.011 



0.005 
0.004 
0.004 
0.003 
0.004 
0.004 
0.003 
0.004 
0.005 
0.002 
0.004 
0.003 
0.003 
0.003 
0.002 
0.002 
0.004 
0.002 
0.003 
0.005 
0.004 
0.004 
0.003 
0.002 
0.003 
0.003 
0.003 
0.003 
0.004 
0.003 
0.003 
0.004 
0.003 
0.002 
0.004 
0.004 
0.003 
0.003 
0.003 
0.005 
0.004 
0.005 



Table 3. continued. 
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0.546 
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0.237 
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0.315 
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0.243 
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0.209 
0.445 
0.345 
0.358 
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0.383 
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0.335 
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0.415 
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0.216 
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0.466 
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0.170 
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0.143 
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0.176 
0.188 
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0.119 
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0.216 
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0.179 
0.258 
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0.104 
0.193 
0.188 
0.124 
0.114 
0.133 
0.086 
0.111 
0.126 
0.080 
0.188 
0.127 
0.136 
0.126 
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0.089 
0.132 
0.088 
0.107 
0.162 
0.162 
0.128 
0.108 
0.131 
0.166 
0.150 
0.118 
0.168 
0.090 
0.113 
0.160 
0.113 
0.087 
0.130 
0.189 
0.101 



0.079 
0.143 
0.140 
0.092 
0.086 
0.102 
0.064 
0.085 
0.095 
0.057 
0.141 
0.094 
0.103 
0.094 
0.087 
0.066 
0.100 
0.067 
0.078 
0.117 
0.121 
0.099 
0.082 
0.095 
0.124 
0.112 
0.085 
0.124 
0.069 
0.081 
0.118 
0.085 
0.065 
0.095 
0.138 
0.075 



0.059 
0.108 
0.106 
0.066 
0.063 
0.080 
0.048 
0.063 
0.071 
0.042 
0.106 
0.066 
0.073 
0.069 
0.067 
0.050 
0.073 
0.047 
0.058 
0.083 
0.091 
0.075 
0.061 
0.069 
0.093 
0.084 
0.063 
0.091 
0.050 
0.059 
0.086 
0.064 
0.048 
0.070 
0.100 
0.055 



0.043 
0.079 
0.079 
0.046 
0.047 
0.062 
0.033 
0.048 
0.053 
0.029 
0.077 
0.046 
0.056 
0.050 
0.048 
0.036 
0.053 
0.030 
0.038 
0.058 
0.068 
0.054 
0.044 
0.049 
0.069 
0.061 
0.044 
0.066 
0.039 
0.044 
0.062 
0.046 
0.035 
0.051 
0.072 
0.040 



0.025 
0.057 
0.060 
0.032 
0.034 
0.048 
0.022 
0.036 
0.038 
0.020 
0.054 
0.033 
0.039 
0.036 
0.034 
0.023 
0.037 
0.021 
0.026 
0.040 
0.048 
0.039 
0.031 
0.032 
0.052 
0.043 
0.031 
0.047 
0.028 
0.031 
0.043 
0.033 
0.024 
0.035 
0.052 
0.026 



0.016 
0.040 
0.043 
0.021 
0.024 
0.035 
0.015 
0.024 
0.022 
0.014 
0.036 
0.022 
0.026 
0.024 
0.024 
0.013 
0.025 
0.014 
0.018 
0.027 
0.032 
0.025 
0.020 
0.023 
0.037 
0.029 
0.020 
0.033 
0.020 
0.020 
0.030 
0.021 
0.014 
0.023 
0.035 
0.017 



0.010 
0.026 
0.024 
0.013 
0.014 
0.022 
0.009 
0.013 
0.014 
0.009 
0.022 
0.014 
0.017 
0.015 
0.016 
0.008 
0.016 
0.009 
0.011 
0.017 
0.019 
0.015 
0.013 
0.015 
0.024 
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0.012 
0.022 
0.012 
0.012 
0.019 
0.014 
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0.014 
0.022 
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0.006 
0.015 
0.013 
0.007 
0.008 
0.009 
0.005 
0.008 
0.008 
0.005 
0.013 
0.008 
0.010 
0.008 
0.009 
0.005 
0.009 
0.005 
0.006 
0.010 
0.011 
0.009 
0.007 
0.008 
0.014 
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0.007 
0.012 
0.007 
0.007 
0.011 
0.008 
0.005 
0.008 
0.013 
0.006 



0.003 
0.006 
0.006 
0.003 
0.003 
0.004 
0.002 
0.003 
0.003 
0.002 
0.006 
0.003 
0.004 
0.004 
0.004 
0.002 
0.004 
0.002 
0.003 
0.004 
0.005 
0.004 
0.003 
0.004 
0.006 
0.005 
0.003 
0.005 
0.003 
0.003 
0.005 
0.003 
0.002 
0.003 
0.006 
0.003 
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Table 4. Galaxy number counts in /, H, and as function of total magnitude, m, expressed in the AB photometric system. 
Differential counts and Icr-errors, n{m) ± An, are quoted per unit mag per square degree. Cumulative counts, A^(> m), are per square 
degree. 
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